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Abstract. We study one-dimensional, continuum Bernoulli- Anderson mod- 
els with general single-site potentials and prove positivity of the Lyapunov 
exponent away from a discrete set of critical energies. The proof is based on 
Fiirstenberg's Theorem. The set of critical energies is described explicitly 
in terms of the transmission and reflection coefficients for scattering at the 
single-site potential. In examples we discuss the asymptotic behavior of 
generalized eigenf unctions at critical energies. 



in L 2 (R). The single site potential f is assumed to be real, supported in 
[—1/2,1/2], locally in L , and not identical to (in //-sense). The coupling 
constants q n {^), n £ 2, are independent, identically distributed Bernoulli ran- 
dom variables, i.e. they have distribution /i with supp \x = {0,1}. Thus the 
random potential Vw(x) = J2 n Q.n{ui)f(x — n) is in Lj^^JM) for all u>, which 
allows to associate a unique selfadjoint operator with H w that may equivalently 
be defined by form methods or in the sense of Sturm-Liouville theory. 

Our interest in Bernoulli- Anderson models of this type arises from the fact 
that their spectral properties are not as well understood, by mathematically 
rigorous standards, as those of Anderson models with continuous (or absolutely 
continuous) distribution \i. Furthermore, these discrete distributions are phys- 
ically more relevant as they are modeling the charge numbers of nuclei. 
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1. The Main Result 



We study Anderson-type random Schrodinger operators 
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One of the central objects in the study of one dimensional random operators 
is the Lyapunov exponent. To define it for our model, let gE(n,u) be the 
transfer matrix of 

-u" + V u u = Eu (2) 
from n — l/2ton + l/2, i.e. for any solution of (|2]) one has 

( «(n+l/2) \ _ , ,( u(n-l/2) \ 
\u'{n + 1/2) ) - 9E{n ^>\ u'{n- 1/2) )' 

Also, for n G N let [/^(n, o>) = gE{n,u) . . . <?e(1, uj), and define the Lya- 
punov exponent at E by 

7 (E) = lim -E(log||[/ B (n,u;)||), (3) 

where E denotes expectation with respect to to. Existence of j(E) follows from 
the subadditive ergodic theorem, e.g. ||. Since \\Ue\\ > 1, we have j(E) > 0, 
and it can also be seen that one gets the same ^y(E) if transfer matrices on the 
negative half line are used analogously in (|3|). 

For many applications it is crucial to know at which energies one has that 
j(E) > 0, which corresponds to exponential growth (or decay) of the solutions 
to (0). Our main goal here is to show that under the above assumptions this 
holds for all but a discrete set M of energies E, and to explicitly describe M 
in terms of the transmission and reflection coefficients for scattering at the 
single site potential /. To define them, let k G C \ {0} and u + (x,k) be the 
Jost solution of 



— u" + fu = k\ (4) 

i.e. the solution satisfying 

u ( xk )-S ^ for x< -1/2, 

+ l ' ' \ a{k)e ikx + b{k)e~ ikx for x> 1/2. { ) 

Since e tkx and e~ lkx are linearly independent solutions of — u" = k 2 u, this 
defines a(k) and b(k) uniquely. They are related to the transmission and 
reflection coefficients used in physics by a = t~ l and b = rt^ 1 . In particular, 
vanishing of b is equivalent to vanishing of r. For real k we use constancy of 
the Wronskian to get 
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|a(*)| 2 -|&(AOI 2 = l, (6) 

corresponding to the familiar |r| 2 + |t| 2 = l. 
It can also be seen from (|5]) that 

/ a(k) \ J_ _ 4fc/2 ( ik -1 \( u+(l/2,k) \ 

V b(k) J 2%k v %keik ~ eik J V w+(V 2 ^) / ' 

Since u + is a solution of (f|), a linear differential equation which is analytic in 
k, and satisfies the analytic initial condition l/2),u' + (— 1/2)) = 

we see that w + (l/2, fc) and u' + (l/2, k) are entire in fc. Thus a(fc) and 6(fc) are 
analytic in C \ {0} with a possible pole at k = 0. Neither a(fc) nor vanish 
identically For a(fc) this is a trivial consequence of (|6|). For b(k) it follows 
from the fact that a compactly supported / ^ can not be reflectionless (i.e. a 
soliton), as follows from inverse scattering theory, e.g. 0, see also [|14|| . Thus 
the roots of a(k) and b(k) can not accumulate, neither at (the pole) nor away 
from 0. 

We are now ready to state our main result: 

Theorem 1. Let E G K\ M, where the set of critical energies M is given 

by 

M := :nEZ}U{E = k 2 :k>0,b(k) = 0} (7) 

U{— a 2 < : a(ia)a(— ia)b{ia)b{— ia) = 0}, (8) 

then i(E) > 0. 

Note that the critical set M is the union of the discrete set {(n7r/2) 2 : n e 
Z} and all real numbers in {k 2 : k G C \ {0}, a(fc) = or = 0} (this uses 
that for real k, a{k) ^ and b(—k) = b(k)). Therefore M is discrete. 

Our main tool for proving Theorem |I] in Section ^| will be Fiirstenberg's 
Theorem, which has been used extensively in proofs of positivity of 7 for dis- 
crete one dimensional random operators, see ||12|, Sec.l4A] for a summary. For 
continuum models, the use of Fiirstenberg's Theorem until recently seems to 
have been restricted to the special cases / = X[-i/2,i/2] an d / = do, a 5-point- 
interaction [jl], §]. The first to have used Fiirstenberg's Theorem systematically 
for continuum Anderson models with general classes of single sites / has been 
Kostrykin and Schrader |9[]. While they state their general results for abso- 
lutely continuous distribution /1, they point out that most of their ideas can 
also be applied to discrete distributions. Our work here can be seen as an 
implementation of this fact with, as we feel, a minimal amount of technical 
effort. 
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A nice feature of using Fiirstenberg's Theorem in the proof of Theorem |T| 
is that once we know that the Theorem holds under the assumption that 
supp fi = {0,1}, then it is immediately clear that the Theorem extends to 
the case where {0, 1} C supp /i. In fact, for larger support the critical set M 



should be smaller. Kostrykin and Schrader note in [10] that M should be 
empty if supp /i has at least one non-isolated point, but there doesn't seem to 
be a proof of this yet. While Theorem [I] does not state that 7(E) = for all 
E G M (which in fact is not generally true) , we will demonstrate in Section |3] 
that for Bernoulli- Anderson models one can indeed find many critical energies 
with j(E) = 0. In the discrete case this has already been observed for the 
so-called dimer model, see 0. In particular, Section ^] discusses the example 
/ = Axr-i/2,1/2]) A G R a constant, which leads to two different types of critical 
energies. They can be classified by the asymptotics of solutions of @. In some 
cases they are of plane wave type, in particular bounded, while other critical 
energies lead to solutions which grow like exp y/x, due to a connection with 
random walks. 

In H we use a result much like Theorem [l] to prove exponential and dynam- 
ical localization for continuum Bernoulli-Anderson models, where in addition 
we build on ideas which were developed for discrete models in |4| . Discreteness 
of the critical set M is crucial to this approach. In particular, it shows that 
the result found in Theorem [1] is superior to the well known results of Kotani 
theory which prove that j(E) > for almost every E in a class of models 
containing ours. The methods of |6[ also allow to extend Theorem [1] to the 
case where the support of /1 is {a, b}, where a 7^ b are arbitrary real numbers. 
This requires considerably more technical effort, due to our use of scattering 
theory at a periodic background, and thereby somewhat obscures the rather 
simple ideas which we present here in Section |2|. 

2. Positivity of the Lyapunov Exponent 

To prove positivity of 7, we need to understand properties of the transfer 
matrices. The transfer matrix from -1/2 to 1/2 of 

— u" + fu = Eu, (9) 

is the matrix, g(E), for which 

u'(l/2) J ~ 9{tj) \ u'(-l/2) J (1U) 

for any solution u of (H). By go{E) we denote the corresponding transfer matrix 
of — u" = Eu. Set G(E) to be the closed subgroup of SX(2, R) generated by 
go(E) and g(E). Let P(R 2 ) be the projective space, i.e. the set of the directions 
in R 2 and v be the direction of v G R 2 \ {0}. Note that SL(2, R) acts on P(R 2 ) 
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by gv = gv. We say that G C SL{2, R) is strongly irreducible if and only if 
there is no finite G-invariant set in P(R 2 ). 

In order to prove Theorem [I] we will consider energies E > and E < 
separately. Note that E = is contained in M, so we don't need to consider 
it. 

The general approach is as follows: We will first prove that G(E) is not 
compact by showing that a sequence of elements has unbounded norm. This 
argument will be valid for all E G M \ M', where 

M' = {(mr) 2 : n G Z} U {E = k 2 : k > 0, b(k) = 0} 
U {—a 2 < : a(ia)a(— ia)b(ia)b(— ia) = 0}. 

Once non-compact, then the group G is known to be strongly irreducible if 
and only if for each v G P(M 2 ), 

#{gv:geG}>3, (11) 

see 0. In order to prove that this condition is also satisfied we will in addition 
have to exclude A;'s which are odd multiples of tt/2. We then use Fiirstenberg's 
Theorem which states, in our context, that if the group G(E) is not compact 
and strongly irreducible then j(E) > 0, see also @. 

We start with positive energies, i.e. E = k 2 , k > 0, and express the transfer 
matrices over [—1/2, 1/2] in terms of the Jost solutions (||]), at ±k. We have 

f u N (l/2,k) u D (l/2,k) \ 
— ^ u ' N (l/2, k) u' D (l/2,k) J 

where and ur> are the solutions of (f|) satisfying u A r(— 1/2, k) = u' D (—l/2, k) 
= 1, u' N (— 1/2, k) = Ud{— 1/2, k) = 0. Writing and ud as linear combina- 
tions of u + (x, k) and m_(x, k) = u + (x, k) and setting z±(k) = a(k)e lk ±b(k)e~ lk , 
we see that 

n (^\-( M z +( k )} lM z +( k ))\ n9 A 

9[K ' - \ -klm[z-(k)] Re[z-(k)} J ' 1 ] 

Clearly, 

n 2\ ( cos k \ sin k \ M oX 

9 ^={-ksink "cask )■ ^ 

Note for k > 0, 

ODU^X^H^)- (14) 

Thus proving non-compactness for G(k 2 ) is equivalent to proving non-compact- 
ness for the group G{k 2 ) conjugate to G{k 2 ) via (|i~4D. 
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If b(k) = we have that \a(k)\ = 1, and hence, we may write a(k) = e 1 ^. 
By (|12"D, one has then that 

w, 2 \ _ ( cos(k + 4>) sin(ft + 4>) \ 
9[ ' ~ \ - sin(A; + 0) cos(A; + 0) )> 

which is a rotation. The same holds for go(k 2 ) by (|13|). Thus G(k 2 ) is a group 
of rotations and thereby compact. This not only excludes an application of 
Fiirstenberg, but shows ^(k 2 ) = as a direct consequence of the definition (|3]). 

Now assume b{k) ^ 0. Set a{k) = Ae ia , b{k) = Be ip , and see that (§) 
implies A 2 - B 2 = 1. It follows then that A > B > 0, A + B > 1, and 
A — B < 1. With this notation, one may recalculate 

= A ( COS ^ +B ( COs( 'f), "I , 

v y V Sin((y9j COS((y9j y V — sm(-0) — cos(^) y 

where we have set := —A; — a and ip := k — (3 . 

We wish to show that a sequence of elements in G(k 2 ) has unbounded 
norm. To do so, consider an arbitrary unit vector 

== ( 2$ ) • 

The relation 

g(k 2 )v(6) = Av(<p + 6) + Bv(-ip - 6) 

shows that the image of the unit circle under g is an ellipse, with the choices 
7] := —\(ip + ip) and rf := | — + -0) identifying the semi-major and semi- 
minor axes as follows: 

g(k 2 )v( V ) = (A + B)v(V=p) and g(k 2 )v(r,') = (A - B)v{\ + . 

Defining R{6) := ||^(A; 2 )f (6*)||, a short calculation, using A 2 — B 2 = 1, yields 

R(6) - 1 = 2B [B + A cos(2fl + <p + $)]. 

As a consequence, we see that R(0) — 1 is 7r-periodic and has exactly two roots 
in [0, 7i"). In particular, R{6) = 1 if and only if cos(26> + cp + ip) = — ~, which 
shows that the distance between the zeros of R(0) — 1 is not equal to |: recall 
B ^ (and B ^ A). Similarly, R{6) > 1 if and only if cos (26 + <p + if>) > -f 
and hence |{0 G [0, tt) : > 1}| > f • As a result, there exists a compact 
interval i\" (not necessarily in [0, 7r)) with \K\ > | and ||g(fc 2 )t> (0) || > c > 1 
for all e K. 

Now, applying g once to v(r]), r\ as above, produces a new vector with 
norm greater than one, but the direction, initially rj, is possibly altered. A 
vector with this new direction may not increase in norm by directly applying 
g again. However, as long as k is not an integer multiple of ir, then finitely 
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many applications of go{k 2 ), which is rotation by k, produces a vector with 
direction in K. Once in K, a direct application of g does increase the norm 
size uniformly by c > 1, as indicated above. In this manner, we can produce 
a sequence of elements in G(k 2 ) with unbounded norm. Thus G(k 2 ) is non- 
compact if b(k) 7^ and k is not an integer multiple of 7r. 

To complete the proof of Theorem [l] for positive energies, we note that if 
in addition k is not an integer multiple of tt/2, then the free transfer matrix 
produces three distinct elements in projective space, i.e. ([D]) is satisfied and 
Fiirstenberg's Theorem implies ^(k 2 ) > 0. 

We finally turn to energies E < 0. In this case G(E) is non-compact for 
every E < 0, since the free transfer matrix 

a(E)=( cosh ( a ) a sinh («) "\ 
^ ' \ asinh(a) cosh(a) J 

has unbounded powers, where a = \/\E\. It remains to check ([TT|), which we 
have to do for E G (— oo,0) \ {—a 2 : a(ia)a(—ia)b(ia)b(—ia) = 0}. 

Let E = —a 2 such that a(±ia) ^ and b(±ia) ^ 0. By (||), this is 
equivalent to 

g(E)v ± £ {v +7 v-}, (15) 

where v± are the directions of v± := (1, ±a)*. If tJ ^ then #{g (E) n v : 

n G Z} = oo. If, on the other hand, v G {v + ,U-}, then we use (p"5|) to con- 
clude that an initial application of g{E) followed by iteration of go{E) gives 
an infinite orbit. This shows ( |TTD and completes the proof of Theorem |1|. 

3. Critical energies 

In this section we discuss the appearance of critical energies in concrete 
examples. We thereby illustrate the following points: (i) critical energies with 
vanishing Lyapunov exponent do indeed exist, (ii) the structure and "size" of 
the set of energies with b(k) = depends strongly on the concrete example, 
and (iii) at critical energies at least two different types of non-exponential 
asymptotics of solutions can be observed. 

3.1. Zero reflection at degenerate gaps. Our first general observation 
is that if H^ er is the Schrodinger operator with 1-periodic potential V satisfying 
V(x) = f(x), —1/2 < x < 1/2 and H^ CT has a degenerate gap at energy 
E = k 2 , i.e. either all solutions of (f|) satisfy periodic boundary conditions 
on [—1/2,1/2] or all solutions satisfy anti-periodic boundary conditions on 
[-1/2, 1/2], then g(k 2 ) = I or g(k 2 ) = -I. By (p)) this implies (a(k),b(k)) = 
(1,0) or (—1,0), respectively. Conversely, (a(k),b(k)) = (±1,0) necessarily 
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requires a degenerate gap. Thus, the possibility of degenerate gaps in periodic 
potentials is one mechanism which leads to reflectionless energies. Whether 
all reflectionless energies arise in this way is equivalent to deciding if b(k) = 
necessarily leads to a(k) = ±1. In Example 2 below we will see that this is 
not the case. 



3.2. Example 1. Consider first / := Axr_i/2,i/2] with A G 1 For E 
k 2 > max{0, A} we get with a := y/E — A 



/ 1 \ ik , ■ 2k 2 A . 

a(k)e =cosa + t sin a 

V ' 2ka 

b(k) = — — sin a 
v ; 2ka 

— ( cos ^ I s ^ n ^ 
' \ —A; sin A; cos/c 



I cos a - sin a 

- i — ' a 

-a sin a cos a 



We can now distinguish between the following types of critical energies: 

Type la: k 2 such that b(k) = 0: This happens for E = k 2 = (rnr) 2 + \, n£N. 
In this case we have g(k 2 ) = ±1, and it easily follows that the group (go, g) is 
bounded. This means that solutions to (0) are bounded and that ^{k 2 ) = 0. 

Type lb: k = nn, n E N: This means Qoik 2 ) = ±1, and again (g , g) bounded, 
7(fc 2 ) = 0. This type is "dual" to Type 1, since k 2 = n 2 ir 2 can be viewed as 
the energies where —d 2 /dx 2 + A — Ax[_i/ 2 ,i/2] is reflectionless with respect to 
-d 2 /dx 2 + A. 

Type 2: E — k 2 , where k = (2n — 1)tt/2, nGN. This leads to a critical energy 
with a rather different asymptotic behavior of solutions if also a = (2m— l)7r/2, 
m G N, that is, for specific values of A = k 2 — a 2 = 7f 2 (n — m)(n + m — 1). In 
this case we have 
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Let h n (uj) := g2 n +i(u)g2n(w), where g n {uS) is the transfer matrix of from 
n-1/2 ton + 1/2. Then 



h n (u) 




with probability p 2 + q 2 , 
with probability pq, 



with probability qp. 

If /i nj i and /i ni 2 are the diagonal entries of h n , then the sums of log \h n j\, 
i = 1,2, give a symmetric random walk (observing that (p 2 + g 2 )logl + 
pq\og(k/a) + qp\og(a/k) = 0). Thus it follows that for every e > and 
almost surely 



lim sup n~ 1/2 ~ £ log || JJ h k (u) || = 



k=l 



and 



limsupn l / 2+e log || J^J hk{oS) \\ = oo, 



fe=i 



see e.g. ||13|1 . The interpretation of this is that the hull of the (oscillatory) 
solutions of (H) asymptotically grows like expc^/x or decays like exp(— Cy/x) 
for some c > 0, where a generalization of the Ruelle-Osceledec Theorem from 



11, Section 81 is used. 



3.3. Example 2. Consider / := A(xr_ 1 / 2j o] — X(o,i/2l) with A > 0. For 
E = k 2 > A, we may again explicitly calculate the transfer matrix 

and see that, 

9n = cos(^/2)cos( Q+ /2)--sin(a_/2)sin(a + /2), 

gi2 = — sin(a_/2) cos(« + /2) H cos(a_/2) sin(a + /2), 

92i = — a+ cos(a_/2) sin(a; + /2) — ct_ sin(a_/2) cos(a + /2), 
g22 = sin(a_/2) sin(a + /2) + cos(a_/2) cos(a + /2), 
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for a± := y/E ± A. From this and the boundary conditions of the Jost solution, 
i.e. 

/ e -ik/2 \ / a (fy e ik/2 + 6 (fc) e -i*/2 

9 y ike -ik/2 J - y ik a (k)e tk l 2 - ikb{k)e~ ik ' 2 

one can show that 

b{k) = if and only if sin(a; + /2) = = sin(a_/2). 

Therefore, if we fix some A > 0, then any E = k 2 > A for which b(k) = 
must satisfy 

= rnr, and = mir, 

2 ' 2 

for some n, m G N. Thus A must have the form 

A = 2n 2 (n 2 - m 2 ), with n, m E N, (16) 

and the corresponding reflectionless £"s are then 

£ = 27r 2 (n 2 + m 2 ). (17) 

From this we conclude that for most values of A (e.g., Lebesgue almost 
every A), no reflectionless energy exists above A. On the other hand, every A 
satisfying fllTf ) has only finitely many reflectionless energies. In fact, if N :— 
2^2 £ N, then N = n 2 — m 2 for only finitely many pairs (n,m) G N 2 . But 
this finite number can be arbitrarily large, depending on A: To find at least j 
exceptional energies, choose 



2tt 2 

and see that for < £ < j one may choose 



Nj = V +1 (2 ] - 1 + 1) 



n e := (2 e (2 J - 1 + 1) + 2^ e - 1 ) 



and 



(2\2 j ~ l + 1) - 2 j - i ~ 1 ) 



satisfying nj — mj = Nj. Proof: Take a := 2 e (2 j 1 + 1) and b := 2 j 1 1 and 
see that (a + b) 2 — (a — b) 2 = Aab. Therefore, 

Aab = 2 2 ■ 2 e {2 j - 1 + 1) ■ 2 j -^ = 2 j+1 (2 j ~ 1 + 1) = Nj. 

Recall from our remarks in ^TT| that in some cases, the condition b(k) = 
implies that a(k) = ±1. To see that this is not always the case, we calculate 
a(k) in this example for such values of k. For A as in (|16D and E as in (|17|), one 
can show that a(k) = e ~ 2m ^ n2 + m2 if n + m is even and a(k) = — e -Wn 2 +m 2 ^ 
n + m is odd. 
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